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tions, and for the second, the limits for small and for large mass are analytically 
obtained too. Coincidence with previously known, partial result serves as a check 
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1 Introduction 



The Casimir effect |IJ is a very fundamental feature common to all quantum field theories, 
which arises in particular, as is well known, when there is a departure from the topology 
of the ordinary flat spacetime towards non-trivial topologies. It has been studied inten- 
sively in the last years due to its importance in elementary particle physics, cosmology 
and condensed matter physics — see Ref. |2j for a recent review on the theoretical and 
experimental aspects of this effect, and also Refs. ||. 

Whenever we deal with a confined charged quantum field, it is a natural question to 
ask for the influence of external fields on their zero-point oscillations and to investigate 
its consequences on the Casimir effect. The influence of external fields on zero-point 
oscillations of unconfined charged bosonic and fermionic fields and the construction of 
the corresponding effective field theories (at the one-loop level) is an honorable subject 
and it is linked to the pioneering works of the thirties by Heisenberg, Kockel and Euler, 
Heisenberg and Euler, and Weisskopf flU and of the fifties by Schwinger ||. More recently, 
the influence of a uniform magnetic field was extensively treated by Dittrich and Reuter 
in the monograph ||, and zeta function techniques were employed to construct effective 
Lagrangians for scalars and Dirac fields for d — 2, 3 and higher dimensions in constant 
background fields (see Refs. [0 and references therein). Ambjorn and Wolfram considered 
the influence of an electrical field on the vacuum fluctuations of a charged scalar field [§]. 
Outside the context of pure QED, the influence on the vacuum energy density of a real 
scalar field due to an arbitrary background of an also real scalar field interpreted as 
a substitute for hard boundary conditions was considered in Refs. ||, [L(| PH . In the 
gravitational case, Elizalde and Romeo investigated the issue of a neutral scalar field in 



the presence of a static external gravitational field (Tj|. A general and valuable discussion 
of these aspects of quantum field theory can be found in Ref. |13|] too. 

The influence of external fields on zero-point oscillations of quantum fields confined 
by hard boundary conditions, however, has been investigated in some particular cases 
only. For example, the influence of a uniform magnetic field on the Casimir effect was 



investigated in Ref. |14| in the cases of a massive fermion field and of a charged scalar 
field both submitted to anti-periodic boundary conditions, and in Ref. |15| in the case 
of a massive charged scalar field submitted to Dirichlet boundary conditions. In the 
cases considered by those authors it was formally shown that the fermionic Casimir effect 
is enhanced by the applied magnetic field while the bosonic one is inhibited. In Refs. 
[HI , 15fl , Schwinger's proper-time method || was employed and explicit analytical results 
obtained for particular regimes of a conveniently defined dimensionless parameter \x and 
magnitude of the external field. In the weak magnetic field regime, some aspects of the 
material properties of confined charged fields were investigated in Refs. flfifl . 

In this work we wish to resume the investigation of the influence of an external mag- 
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netic field on the zero-point oscillations of charged confined quantum fields and consider 
the case of a massive fermion field under a uniform magnetic field and constrained by 
boundary conditions of the MIT type, which states that the fermionic current through an 



hypothetical confining surface must be zero ||17|| . This is not an academic question since, 
for instance, in the bag model of hadrons [|l7], we can expect the zero-point oscillations 



of the quark fields to be influenced by the strong electric and magnetic fields that perme- 
ate the interior of the hadronic bag (see for example and references therein) . We wish 
also to reconsider the case of a massive charged scalar field under a uniform magnetic field 
by investigating regimes not considered in Ref. JT5|. In order to evaluate the relevant 



Casimir energies it is convenient, especially in the case of the massive fermion field with 
MIT boundary conditions, to follow the spirit of the representation of spectral sums as 
contour integrals pO|j , which allows for the incorporation of the boundary conditions in 
a smooth way. Here we will employ a relatively simple variant of this general procedure. 
This variant is described in and it will be here applied without further explanation 
to the cases at hand (the interested reader should consult that paper for details). The 
outline of the present work is as follows. In Section 2 we apply our calculational tools to 
the charged scalar field case. In Section 3, we consider the case of the charged fermion 
field. For both cases, detailed numerical analysis and some special analytical limits are 
provided, together with an explanation of the results obtained and specific comparisons 
with other results. Section 4 is devoted to final remarks. Throughout the paper we employ 
natural units (h = c = 1.) 



2 The vacuum energy of a charged scalar field with 
Dirichlet boundary conditions and under the influ- 
ence of a uniform magnetic field 

Let us first briefly consider a charged scalar field under Dirichlet boundary conditions 
imposed on the field on two parallel planes separated by the distance £ whose side L 
satisfies the condition L ^> I. Suppose initially that there is no external field. The 
unregularised Casimir energy is given by [2TH PI 



E Q (£, /i) = a 



(2*0 



d 3 p log 



1 + 



Ko(z) 



(1) 



where K\ {£) and K% (z) are functions constructed from the boundary conditions as de- 
scribed in Ref. pTJ. The dimensionless parameter a takes into account the internal 



1 Notice that in the absence of external fields we could have started from Eq. (21) in Ref. |2l[], but for 
our purposes here it is more convenient to start from Eq. (18). 
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degrees of freedom of the quantum field. For Dirichlet boundary conditions we begin by 
writing 

F (z) = sin z. (2) 

Since z = is a root of F(z) we divide this function by z, thus removing z = from 
the set of roots without introducing a new singularity: This is equivalent to removing the 
zero mode. Define 

G (z) = !El (3) 

The function K (z) = K\ (z) + K 2 (z), with K\ (z) = K 2 (—z), is then obtained by per- 
forming the substitution z — > iz, that is 

K{z):=G{iz) = ^—, (4) 

where z is the function 



z = z(p lt p 2 ,p 3 ) = i^Jpl+pl+pj + m 2 . (5) 
Since K\ (z) = —e~ z /2z and K 2 (z) = e z /2z we have 

E (£, /i) = [ log [l - e-H . (6) 



2 J (27TY 

To shorten these initial steps let us set m = 0. Then after expanding the log we obtain, 
simply, 

where we have defined x := p^i. The integral can be evaluated with the help of the Mellin 
transform 

A"T(s) = / dtt s - l e~ At . (8) 



and after, simple manipulations, we obtain the well known result 

^ /.x aL 2 7i 2 

For a = 2 we get the result valid for photons. A more complex example of the usefulness 
of Eq. ([I]) is provided in Ref. . 

Let us consider now the same type of field and boundary conditions in the presence 
of an external magnetic field which we will suppose to be uniform and perpendicular to 
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the two Dirichlet planes. We will assume also that eB points towards the positive OX 3 
direction. Equation ([!]) reads now 



E (£, /i, eB) = a 



feB\ tf™ y°° dp 



\2irJ 2 



E 

n=0 



2tt 



log 



1-e 



-22 



(10) 



where we have taken into account that for a charged spin zero boson the Landau levels 
are given by 

p\ + pi = e B (2n + 1) , n = 0, 1, 2, 3, ... . (11) 

The factor eB/2n is the degeneracy factor and for a charged scalar field a = 2. Hence, 
the function z reads now 



z = z (p 3 , n) := ^Fp 2 . + eS£ 2 (2n + 1) + /i 2 , 
where \i := £m. It is convenient to define 



and write 



where 



M n := \J(2n+ l)eB£ 2 + /j, 2 

eBT 2 00 
E (£,^eB) = —J2ln(M n ) 



n=0 



In{M n 



dx log 



1 _ e -2 v /x2+M2 



(12) 

(13) 
(14) 
(15) 



and we have set a; := p 3 £. Expanding the log and introducing the variable u, defined by 



uj := Jx 2 + M 2 , we end up with 



eBL 



2 oo oo 



E (£, /i, eB) = —r £ E r^n (M„ 



n=U fe=l 



where 



hn (M n 



(16) 
(17) 



duo uj (u + M n )~ 1/2 (w - M n )~ 1/2 e- 2fc ". 
In order to evaluate this integral, we first introduce an auxiliary integral defined by 

I kn (M n , A) := [°° duo {uj + M n y 1/2 (uj - M n y 1/2 e- 2k "\ (18) 



Mn 



To evaluate the auxiliary integral we make use of (c./. formula 3.384.3 in [|22 



cla; (# + /?) (x — fj,i) H e^ 2 = — exp 



i-'i 



•f*2 



xT (2p) W u _ p , u+p i + Pto) , (19) 



which holds for fix > 0, |Arg (/3 + //i)| < 7r, Re /i2 > and Re p > 0. The auxiliary 
integral then reads 



7 fcn (M n , A) = ^4u72 r (V2) Wo,o (4fcM n A) , 



(20) 



where W^a (-2) is the Whittaker function J23j. The Whittaker W ,o (z) function is related 



to the modified Bessel function of the third kind through [22 



1 /2 

^0,0 (4kM n X) = (4fcJ ^ } K (2kM n X) . 



71 



(21) 



Taking the derivative of Eq. (p0[) with respect to A, and setting A = 1 to obtain l kn {M n ), 
we have 

I kn (M n ) = 2M n —K (z = 2kM n X) | A=1 = —M n K x (2kM n ) . (22) 
Hence the vacuum energy reads 



E (£,fi,eB) 



eBL 



2 00 



M n T Ki (2kM n 



Z7F * n=0 k=l K 



or, more explicitly, 



(23) 



E (£,fx,eB) 



eBL 2 

' 27i 2 e 



eB£ 2 + p 2 Y, T Kl ( 2k y eBi2 + F 
k=i k 



eBL 
~2t? 



2 CO 



Y, \/( 2n + !) eB£2 + A* 2 H ( 2k V ( 2n + !) e5£2 + A* 



n=l 



fc=l 



(24) 



This representation is an alternative to the one obtained in ||15fl . In a very strong field, 
such that eB£ 2 ^> /x 2 , the decaying exponential behavior of the Bessel functions of the 
third kind allows us to keep only the term corresponding to k — 1 in the first parcel of 
the vacuum energy and, thus, in this limit we have 



E (£,eB) ^ (eB£ 2 f 4 ^ 



L 2 



7T 



3/2 £3 



(25) 



in agreement with |I5| |. Notice that for zero magnetic field the vacuum energy as given 
by Eqs. (j2~lj) and (|25|) is zero, that is, the zero of the energy is automatically shifted with 
respect to the vacuum energy in the absence of the external field (which serves, therefore, 
as the natural origin of energies). 
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2.1 Numerical analysis for the scalar case 

For arbitrary values of the parameter \i and of the scaled field eB£ 2 , a numerical eval- 
uation of Eq. (|23|) is still possible. The graphs of Figs, la, lb, 2, and 3 exhibit some 
representative examples. 




Figure 1: Plot of the energy in terms of the dimensionless quantity b(= eBl 2 ) corresponding to 
the magnetic field, for a fixed value of the dimensionless mass /i (here fj, = 1). Fig. la shows in 
detail the formation of a smooth minimum and its precise value. 



In the first two graphs (Fig. 1) we plot a dimensionless version of Eq. ( f23| ) as a 
function of the scaled magnetic field b := eB£ 2 , for fi = 1 and a convenient range of b that 
clearly shows the behavior of the energy (Fig. lb). Notice that initially the external field 
decreases the value of the vacuum energy up to a certain value of b for which the vacuum 
energy attains a minimum (which is most clearly depicted in Fig. la). After reaching 
this point, the energy grows as remarked in Ref. [15|]. This situation is the opposite of 
what happens with a fermionic field, in which case the energy decreases linearly with b, 
for large values of b (as will be clearly seen in the next section). 

Fig. 2 shows the behavior of the energy in terms of the mass \i for a fixed value of 
the magnetic field b (chosen here as b = 1, but the shape of the curve is very similar for a 
wide range of values of b). It starts at a non-zero value, for zero magnetic field. Also this 
is quite different from the behavior in the case of a fermionic field. 

Finally, in Fig. 3 the energy is depicted vs the variable t = b/fi 2 , that measures the 
magnitude of the dimensionless magnetic field b in terms of the dimensionless mass /i 2 . 
Notice that, although for \i = 1 we obtain the same curve as in Fig. 1, for any other value 
of \i the precise shape of Fig. 3 is valuable, as a way of representing the energy in terms 
of the unique quantity t (the magnetic field in units of mass squared) . 

Families of curves corresponding to different values of one of the variables, and also 
two-dimensional graphs are easy to obtain, in a reasonable amount of time, from our 



7 



E 




mu 



Figure 2: Plot of the energy in terms of the dimensionless mass fx, for a fixed value of the 
dimensionless magnetic field b (here 6 = 1). For a wide range of values of b one gets a curve with 
a similar shape. 



formulas in this section. 



3 Confined fermion field in a uniform magnetic field 

As in the case of the charged bosons, the uniform magnetic field is here perpendicular 
to the (hypothetical) parallel MIT constraining surfaces. The distance between them is 
£ and eB is again supposed to point towards the positive OX 3 direction. The starting 
expression reads 

1 ( eBL?\ 00 

E Q {!,n,eB) = -2x-x[—— £ £ I na , (26) 

1 \ A7T / n=0 ae{-l,l} 

where the multiplicative factor 2 takes into account the particle and antiparticle states, 
and we have defined 

T r d P 3 1 

La ■= / — lo; 

J-oo Z7T 

where again (with \i := £m) 

z = z(q, n, a) := y^ 2 p| + (2n + 1 - a) eB£ 2 + fj, 2 n = 0, 1, 2, 3... . (28) 



Z + pL 



(27) 
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Figure 3: Plot of the energy in terms of the dimensionless variable t = b/ fi 2 , that measures the 
magnitude of the dimensionless magnetic field b in terms of the dimensionless mass fi 2 . It is clear 
that, at fi = 1 we recover the same shape as in Fig. 1. 




Expanding the log as before we obtain 
log 



1 - Z -^e~>* 



z + fl 



°° (-1 

E 



k=l 



k 



x [z — fi] e 



k -2kz 



and performing the sum over a, we arrive at 



eBL 



2 oo ' oo 



£ (^, e 5) = -2x— £ £ 



where 

/pfc (M„) : = 
with x := p 3 £ and 



c/.r 



(x 2 + M p 2 ) 1/2 + /i 



p=-l k=l 



x 2 + M; 



1/2 



e -2fc(x 2 +M p 2 



,V2 



M 2 := 2 (p + 1) e5£ 2 + fi 2 , p= -1, 0, 1, 2, 3, 



(29) 



(30) 



(31) 



(32) 



The prime in Eq. (|30|) means that the term corresponding to p = —1 must be multiplied 

1/2 



by the factor 1/2. Let us define, as before, a new variable u through u = yx 2 + M 2 
Then 

I pk (M p ) = [°° du u (u + M p y 1/2 (uo - M p )~ 1/2 (u + fi)~ k {uo - fif e~ 2ku . (33) 

J Mr, 



The integral defined by Eq. (|33|) is non-trivial. Here we will solve it analytically in two 
limits, to wit: the large- and the small-/x limits. A numerical integration of Eq. (|33|) , 
however, is feasible to conveniently complete the analysis. We will show the results later. 



3.1 The limit 

If we set /i ~ we have to solve a much simpler integral for l pk (M p ). In fact, 

f°° uje~ 2ku} 
I pk (M p )= du (34) 
Jm p ^Juj 2 - Mj 

This integral can be evaluated with the help of (cf formula 3.365.2 of fl22|) 

/ = aKi jab) , a > 0, Re a > 0. (35) 

J a vi 2 — a 2 

For a = this integral reads 

oo 1 
bx L 







e~ bx = v (36) 

as can be easily seen by recalling that, for z — ► 0, K\ (z) wl/z, There are three types of 
contributions to the vacuum energy: corresponding to p — — 1, p — 0, and p = 1, 2, 3, 
respectively. They read 

eBL 2 ~ (-l) fe+1 r e£L 2 ~ (-l) fc+1 r 

e5L 2 ~ ~ (-l) fc+1 



7T 



2/ 



E E , (37) 



p=l k=l 

where the integrals in the partial sums are given by 



I-iu = Y k > ( 38 ) 



I ok = V2eBi 2 Kt (2kV2eBe 2 ) , (39) 

L pk = yj2 (p+^eBi 2 ^ [ik^jl (p + 1) eB£ 2 ) . (40) 

It follows then that the first sum can be exactly evaluated by using the eta function, 

Vr ( z ) = (1 — 2 1 ~ s ) ( R (s). In fact, we have 

oo / _2 
fe=l 
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Therefore, the Casimir energy in this limit reads 



2 oo / 



p rjr2 pT^T 2, 00 ( 1 \ 



48 1 ir 2 £ £-! k 

k=l 



eBL 2 ~ ~ (-l) fe+1 



, , y: E ^^V 2 (P + !) e5 ^ 2 ^ 2 V 2 (P + !) e ^ 2 • ( 42 ) 

71 * p =i fe=i K v y 

From this result we can easily extract the limit eB£ 2 3> 1. Due to the decaying exponential 
behavior of the modified Bessel function, the leading contribution is given by the first term 
in Eq. §2]). Hence 

eBA 

£*,(^«0,eB)«-— , e5£ 2 »l. (43) 
This result has the same sign — and one fourth of the magnitude — of the corresponding 



problem with anti-periodic boundary conditions in this same limit []T4|]. This difference 
can be readily understood if we recall that the anti-periodicity length corresponds to 2£, 
and that the number of allowed modes is twice the number of modes of the MIT case for 



massless fermions. Therefore, the result given by Eq. (|43|) , in that particular limit, is 



compatible with the one given in Ref. []T4 . 



3.2 The limit ji » 1 

In this limit only the term corresponding to p = —1, M_i = fi, contributes. Therefore, 
we have to calculate 

I_ lk fa) = T duu{u + iiY k - l/2 (to - fx) h - 1/2 e~ 2kuJ (44) 

The trick to evaluate this integral we have learned it above. One should write the auxiliary 
integral 

I_ lk fa, A) := I" duo {uo + fi)- k - 1/2 (u - fi) k ~ 1/2 e~ 2ku)X (45) 



and evaluate it with the help of formula given by Eq. (|19|) . The result is 

I-ik A) = T ^^ W. m (4/ifcA) . (46) 
(4/x/cA) 

For large values of the argument, the Whittaker function behaves as 

W ull (z)^z^e- z / 2 J (47) 



11 



hence for // ^> 1 we can write 



J_ lfe (/x, A) « T (fc + 1/2) 



-2/ikX 



(4/xfcA) 



fc+1/2 ' 



and since 



we finally have 



J - lfc(/i ' A) = -2l^ 



r (fc+1/2) 



J_lfc (/i, A) 



r (fc + 1/2) 

22fc+2 fck+3/2 



(4^A) fc+1/2 
(A; + 1/2) 



fc+1/2 



+ 



2fc 

fc-1/2 



A=l 



3 -2^fc 



(48) 



(49) 



(50) 



The vacuum energy in this regime is given bay 



pRT 2 00 

^(i jA4 >i, e s)«-i- 57 . x: 



\k+l 



8vr 2 £ ^ 2 2k k k * 5 / 2 



r (fc + l/2) 



(fc + l/2) 2fc 



fc+1/2 



fe-l/2 



-2fik 



, (51) 



a result that is new and incorporates the first mass corrections. 
The most relevant one is the fc — 1 term , which yields 



eBL 2 e" 2 " 
"32tt 3 / 2 £ fi 1 / 2 '' 



(52) 



that is sufficiently small. This is not the same result as the one obtained for anti-periodic 
boundary conditions though the damping exponential appears in both cases. The factors 
multiplying the exponential are different, it must be remembered though that for fi 
the MIT and the AP spectra are not comparable. When /i = the MIT and the AP 
spectra for each component of the fermion field differ by a factor of 4. 



3.3 Numerical analysis for arbitrary mass 

For arbitrary values of the mass and the magnetic field or, correspondingly, of their 
dimensionless counterparts, \i and eB£ 2 , the integral defined by (|33|) is very hard to 
evaluate analytically. In fact, even if we would give it in terms of hypergeometric functions 
this would not improve our knowledge of the dependence of the energy on the mass and 



the magnetic field. Fortunately, a numerical analysis based on Eq. (|30|) is possible and 
leads to very precise, easily understandable results. 

To this end, first we rewrite Eq. (|30|) under the form 



£ 3 E (t, ft) tf£ 



3 oo ' oo / i \fc+l 



4EE y -^—J pk {t,ii), (53) 



L 2 7T 2 

n H P =-i k=i 
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where 

/ 9 , x \V2 l fe 

(y 2 + 2(p + l)t + l) -1 

xe -2k fl (y2+2(p+l)t+l) 1/2 ^ 

with y := x/fi and i = eB£ 2 / ji 2 . 
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Figure 4: In the fermionic case, plot of the energy in terms of the dimensionless magnetic field 
6, for a fixed value of the dimensionless mass \x (here fi = 1). In Fig. 4a, detail the formation of 
an inflexion region, and the inflection point; in Fig. 4b, the intermediate region; finally, in Fig. 
4c, asymptotic behavior for large values of b. 

The result of the numerical evaluation of Eq. (|54|) for a sample of values of /i and b 
is shown in Figs. 4 and 5. Fig. 4 shows the dependence of the energy on the magnetic 
field b for different ranges of this quantity, to better show the inflexion region (Fig. 4a, 
low values of /z), the intermediate region (Fig. 4b) and the asymptotic behavior for large 
values of b (Fig. 4c). The intermediate region is most interesting, showing a smooth 



dy 



y 2 + 2 (p + 1) t + 1 



1/2 
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Figure 5: Plot of the energy in terms of the dimensionless mass fx. In Fig. 5a we see the behavior 
for small mass, with the formation of a smooth minimum, which looks much sharper in Fig. 5b, 
where the asymptotic behavior for large fi is clearly established. 

transition from the case of small magnetic field (with a good zero field limit) to the 
linear asymptotic behavior corresponding to a large field. And this situation is common 
for any particular value of the mass ji. The transition is given by a smooth inflexion 
point of nearly horizontal tangent. It looks like a minimum could almost be formed, 
but not quite. In Fig. 5 the dependence of the energy on the dimensionless mass ji is 
depicted, showing the behavior in the region of small mass (Fig. 5a) and the asymptotic 
behavior for large \x (Fig. 5b). Needless to say, the asymptotic behavior corresponds to 
the the analytic expressions obtained before. For the intermediate region we see a smooth 
behavior connecting the two regions with a minimum that is very easily obtained with 
good precision, for any particular value of b. 

The same observations as for the bosonic case can be done here, namely, that the 
graphs shown are only a very small sample of the ones that can be easily obtained from 
our general formulas in this section, corresponding to the fermionic case, for different 
values of the variables, and including two-dimensional plots. We should warn the reader, 
however, that the computation time is drastically increased now. 

4 Final remarks 

In this paper we have considered the influence of an external uniform magnetic field on the 
Casimir energy associated with charged quantum fields confined by hard boundary condi- 
tions. The method employed incorporates, in a relatively simple way, the hard boundary 
conditions of the Dirichlet or MIT type imposed, respectively, on a massive bosonic field 
and a massive fermionic field plus, for both cases, the effect of the magnetic field, and 
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leads to expressions for the vacuum energy especially suited for numerical calculations. 
The analytical and numerical results show that, in the bosonic case and Dirichlet bound- 
ary conditions, the effect of a large, applied magnetic field is to suppress the vacuum 
energy, and in the case of fermionic field confined in a slab-bag with MIT boundary con- 
ditions, the effect of the applied field is to enhance (in absolute magnitude) the vacuum 
energy. For small and for intermediate values of the scaled magnetic field, however, a 
more interesting behavior of the vacuum energy shows up. Since here the calculations 
were performed in the general framework of an effective field theory, the details of the 
interaction between the quantum field and the external magnetic field — represented by 
the generation of the Landau levels and the hard boundary conditions — remains hidden. 
It would be interesting to understand the behavior of the vacuum energy under the cir- 
cumstances considered here in terms of a more fundamental model of the structure of the 
confined quantum vacuum. 
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